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Abstract. Using the density matrix renormalization group (DMRG) method, we study 
the quantum coherence in one-dimensional disordered spin chains and Fermi systems. We 
consider in detail spinless fermions on a ring, and compare the influence of several kinds of 
impurities in a gapless and a dimerized, gapped system. In the translation-invariant system 
a so-called site-impurity, which can be realized by a local potential or a modification of 
one link, increases for repulsive interaction, and decreases for attractive interaction, upon 
renormalization. The weakening of two neighbouring bonds, which is a realization of a so- 
called bond-impurity, on the other hand, is healed for repulsive interaction, but enhanced 
for intermediate attractive interactions. This leads to a strong suppression of the quantum 
coherence measured by the phase sensitivity, but not to localization. Adding a local distortion 
to a dimerized system, we find that even the presence of a single site-impurity increases the 
metallic region found in the dimerized model. For a strong dimerization and a high barrier, 
an additional sharp maximum, is seen in the phase sensitivity as a function of interaction, 
already for systems with about 100 sites. A bond-impurity in the dimerized system also opens 
a small metallic window in the otherwise isolating regime. 
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1 Introduction 

Recent experiments have led to a renewed theoretical interest in disordered spin-Peierls 
systems. For example, the doped one-dimensional Heisenberg system Cui-^Zn^GeOs 
Jy shows two subsequent magnetic transitions, the spin-Peierls transition at T = 14 
K followed by a transition to an antiferromagnetic ground-state at Tn = 5 K, while 
in Cui-^Mg^GeOs, a reentrant spin-Peierls phase for x > x c Q is observed. In fact, 
the general question of the effect of various types of impurities in pure || |J, |j| [| [?J ||] 
or gapped spin systems 10 , |ll) has been studied intensively during the last years. 



Nevertheless, the interplay between interaction, disorder, and periodic distortions is 
still a challenging problem. The model of spinless fermions, which is considered here 
in detail, is equivalent to the anisotropic Heisenberg model. It describes certain as- 
pects of magnetic and electronic systems, and the phase diagram of this "simplest" 
interacting fermion model is surprisingly rich. On the basis of this model, the aim 
of this work is thus to achieve a better understanding of the ground-state properties 
of disordered spin- and interacting Fermi-systems, especially clarifying the role of the 
interaction, which may enhance or decrease the localization due to the random and 
periodic perturbations. In particular, we will introduce and compare the effects of 
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various kinds of impurities. 

In this context, Eggert and Affleck Q pointed out that two "classes" of impurities 
exist, which differ in their effect on the local symmetry of the system. So-called 
"site-impurities" violate the site parity by affecting one site or bond, whereas "bond- 
impurities" , which modify two neighbouring bonds, violate the bond symmetry but 
respect the site parity. Bearing this in mind, we use the name site- or bond-impurity. 
Note that this naming is different from the one used in and || , where an impurity 
on a single site is called site-impurity, and an impurity on a bond is called bond- 
impurity. Another type, "transparent" impurities, which we do not consider here, 
can be constructed using the Bethe-Ansatz Those defects are similar to bond- 
impurities, but include also a coupling between the next-nearest neighbours around 
the impurity. Transparent impurities have also been studied extensively with analytic 
and numerical techniques during the last years, see for example |Q and ||. As a result 
of Eggert and Affleck, a site-impurity is relevant in the sense that it can break up 
a closed ring, whereas a bond-impurity is irrelevant, i. e. the defect is healed at low 
energies. 

Gapped spin systems containing irrelevant impurities, i. e. impurities which do not 
close the gap in the system, are supposed to be equivalent to free Dirac-fermions with 
random mass, a model which has been widely investigated in one dimension Jl2|. 12] 
in the context of doped spin-Peierls or spin-ladder systems and two dimensions 14] 
in the context of the Quantum Hall Effect. However, the transformation to Dirac 
fermions is valid only for special points in the parameter space, which correspond 
to the XY model (free fermions) for the dimerized system, or the isotropic XXX- 
Heisenberg antifcrromagnet for the ladder model. 

The numerical density matrix renormalization group [l5| is a quasi-exact numerical 
method to determine the ground state properties, in particular the ground-state en- 
ergy, of long one-dimensional (non-integrable) systems with reasonable accuracy. The 
bosonization technique ][l6] can be used with advantage to interpret the numerical 
data. 

In the following, we introduce the models and the impurity types studied by us. 
Using the bosonization technique and the Luttinger description, we then identify the 
leading non-linear operators, and classify the impurities according to them. The results 
for single defects are presented in Sec. ||. Concerning the dimerized system, we 
concentrate in our discussion on strong site- and bond-impurities. The results are 
discussed in Sec. ||, followed by a summary in Sec. [|. 



2 Magnetic chains and spinless fermions 



We begin by presenting the generic spin model, and then describe the equivalent 
fermionic model. Then we briefly introduce the Luttinger description, which is useful 
for discussing the relevant operators which lead to insulating behaviour. 



C. Schuster and U. Eckern, Quantum spin chains with defects 



3 



2. 1 The Heisenberg spin chain 

As a starting point for the study of a general disordered spin-Peierls system in one 
dimension, we consider an anisotropic Heisenberg chain, given by the XXZ model, 
with a dimerized interaction: 

N K 
ffspin = -£ Ji{u) (a?<?? +1 + aVaf +1 + Aaf <r? +1 ) + N-u 2 , (1) 

where the dimerization in the Peierls state, u, enters the coupling constant according to 
Jj(ii) = J[l + (—)*«]. For the clean XXZ model, i. e. for u = 0, and for zero total mag- 
netisation, M = = 0' one finds three phases: a ferromagnetic phase for A > 1, 
separated by a first-order transition from a gapless phase for — 1 < A < 1, whose low 
lying excitations are given by those of a Luttinger liquid; and an antiferromagnetic 
phase for A < — 1. The transition from the Luttinger to the antiferromagnetic phase 
is of Berezinskii-Kosterlitz-Thouless type. Adding the dimerization, u, the system 
becomes localized by forming spin-singlets on neighbouring sites for A < 0, i.e. for 
antiferromagnetic coupling. An excitation gap opens due to the usual Peierls mecha- 
nism |]l7f . Leaving aside the question whether a finite u can be stabilized, we note that 
the dimerization is already relevant for A < \/2/2 jl8|, |l9), and the ground-state wave 
function is localized. The interaction-dimerization phase diagram was determined in 

Impurities may be realized in three different ways. First, local magnetic fields 
may be present, resulting for example from magnetic impurities near the chain, which 
couple directly to the er z -component at a specific site. We call this kind of impurity 
a site-impurity. Second, the coupling between the spins may be modified locally by a 
factor (1 — <5), for example by substituting Sr with La or Ge with Si in a spin chain 
compound like La x Sr2_ x Cu03 or CuGei_ x Si x 03. As only one link is concerned this 
is similar to a local field; nevertheless we call it link-impurity in the following. The 
third possibility is to modify both bonds to the left and to the right of the impurity 
site by an equal amount, (1 — 6): This situation, which is called bond-impurity, can be 
realized by doping magnetic impurities with spin 1/2 into a chain, e. g. substituting 
V by Nb. (Nevertheless, in some cases it is found that the free Nb-electron - with the 
spin - moves to the Vanadium site, i. e. a Nb 5+ - and a V 3+ -site is established instead 
of Nb 4+ and V 4+ .) Completing the Hamiltonian ([!]) in this way, we can write 

N N 

H spin = -J2 Mu, 6, b) (a>f +1 + afof +1 + Aaf - £ hrf , (2) 

i=l i=l 

with Ji(u, S, b) = Ji{u)[\ — Si im 5 — (Si im + 5i, rn +i)b], where m is a fixed site. 

2.2 The fermionic model 

The corresponding fermionic model is obtained via the Jordan- Wigner transformation. 
In the result we change the notation, J — > i, JA — > — V/2, and e m — > —2h m , and 
neglect constant energy shifts, like ^ hi and A (2 Nf — A/2), where Nf is the number 
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of fermions: 

-fffcrmion = ~ ^ U (c+ C i+1 + C+^Cj ) + ^ V^Tlj+l + ^ 6 i n i ( 3 ) 
i i i 

v 

-~2 + n m+i) + b( n m-i + 2n m + n m+ i)\ , (4) 

where U(u) = t[l + (-l) l u t - S itm S t - {5^ m + S iym+1 )b t ], and Vi{u) = V[l + {-l) l u v - 
8i t m3v—(fii.m+8i t m+i)bv]- Furthermore, we have assumed that the coupling to (static) 
phonons and to the impurities can be varied independently in both the hopping and 
the interaction, by introducing six parameter (u t , uy, S t , Sy, b t , by) instead of three 
(u, S, b). The particle density is n — Nf/N, and we restrict ourselves to half filling; we 
set t = 1 in some of the formulas below. As shown in pl[ |, u t is the main contribution 
arising from the dimerization, whereas uy modifies the results only quantitatively but 
not qualitatively; thus we assume uy = 0. For simplicity, we likewise neglect (||) to 
avoid the mixing - in the fermionic picture - of bond- and site-impurity. In the case 
of the link-impurity, it was shown by Meden et al. Q that the weakening of one 
bond in the spin model is equivalent to an impurity given by St and 8y, i. e. that 
the contribution (||) is of lesser importance. We assume that the same holds for the 
bond-impurity. The in this way simplified Hamiltonian is then given by 

-fffermion = - ^ t[l + (-l) l U t ~ S iym 6 t ~ (6 i<m + S i>m+1 )b t ] (c 4 + C, i+1 + C+^C, ) 

i 

+ ^2 V[l - 5 ijm 8y - (Si <m + 6 iym+1 )by] + ^ e i n i- ( 5 ) 

i i 

Thus the link- impurity can occur as a pure hopping- impurity, 5 t ^ and Sy = 0, 
or a generic link-impurity, S t ^ and Sy ^ 0. Setting St = 0, we can also realize 
an interaction-impurity by varying Sy. In the following we consider the Hamiltonian 
given by Eq. (||) only. 



2.3 Bosonization and Luttinger description 

For the interpretation of the numerical data, we use the fact that in the gapless phase 
the system can be described by a Luttinger liquid, and that the distortions can be 
considered as a perturbation. In the bosonized form (TJ), the Hamiltonian can be 
written for the "clean" model as follows: 

H hoson = J^- |^ [d x <p{x,t)f +w 5 [7rn y (a;,i)] 2 | , (6) 

where IT^ is the momentum density conjugate to ip(x, t), i. e. it corresponds to 11^ = 
dt<p/(vg) in the (Euclidean) path integral formulation of the theory. It is also related 
to the conjugate phase variable 6 through 11^ = d x 8(x)/-K. The velocity v of the 
bosonic excitations is given by v = [irt sin(2?7)]/(7r — 2r]), and the interaction constant 
is g = 7r/4?7, where r\ parameterizes the interaction according to V — —2tcos(2rj). The 
density (operator) is given by 

d k 

n(x) = no H — I cos [2kpx + 2ip(x)]. (7) 

TTOX 7T 
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We chose the above standard representation because then the order of the scatter- 
ing process is directly seen in the non-linear terms. Non-linear - global and local 
- contributions arise from the dimerization and the interaction as well as from the 
impurities. 

• Dimerization causes a 2fci?~scattering process of the fermions, 

H u oc 2u J dx snx\2y(x)]. (8) 

• While the fermion-fermion scattering with q « and q ~ 2k F is absorbed in the 
Luttinger-parameter g, the 4/cF-scattering leads to 

H v oc V J dx cos [4tp(x) + (4k F x - G)x], (9) 

where G is the reciprocal lattice vector. This term causes the transition to the 
CDW-state for V = 2 (g = 1/2). 

• Backscattering arises from a local potential (say at x = 0). Since the potential 
couples directly to the density, Eq. (Q), we find in lowest order in the impurity 
strength 

H e = e ?io oc e n 2 k F {x = 0) oc e cos[2<p(0)]. (10) 

A strong potential in an infinitely long chain is equivalent to a weak link |^2| : 
The diagonalization of the above Hamiltonian, -fff cr mion + H e , leads, in the case 
of a strong potential between two semi-infinite chains (R: right chain, L: left 
chain), to the contribution 

U t oc C « (0)Cl(0) oc - cos[20(O)]. (11) 

• A change in the hopping is also considered in the weak, St — > 0, and strong, 
St — ► 1, limit. Following the analysis for the (weak) periodic [^l| or random 
potential p3[ , we would assume that the transition from a site-centered potential 
to a bond-centered potential causes only a phase shift of tt/2 in the non-linear 
term, i. e. a shift from a cosine to a sine. The weak link thus should correspond 
to a high barrier, see again i. e. 

H t oc St sin[2p(0)] , S t (12) 
H t oc (1 - St) cos[20(O)] , S t ^l (13) 

should be appropriate descriptions. However, if we consider a modified link with 
H t oc (1 — St)(ciC Q + h. c), we find instead 

H t oc 5tsinfo>(l)Binfo>(0)] (14) 
H t oc (1-J t )cos[0(l)-0(O)]. (15) 



G 
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• The local modification of the interaction, i. e. the interaction-impurity, leads to 
a 2fcp -contribution too, in addition to the 4fc^-term already present in Eq. (||): 

H Sv oc H t (16) 

• The modification of two links, in the limit of b = bt = by — ► 0, does not cause 
2fci?-scattering at half filling. Thus, no non-linear contribution to the Hamilto- 
nian is found. In the limit b — > 1, however, we find the following expression: 

H b = (l-6)(cfc b + cfc_ 1 +h.c) 

oc (1-6) sin{[0(l) + 0(-l)]/2 - 9(0)} sin {[9(1) - 0(-l)]/2} 

xcos[^(0)] sm{[^(l) + ^(-l)]/2} 
oc (1 - 6)sin($/2)sin[9(l)]cos[^(O)]sin[0(l)]. (17) 

In the last step we use 9(x) = 9q + <frx /L + 9(x) , see jl(| , where 9q represents the 
zero- mode contribution, <S>x/L the boundary condition see also the following 
section-, and 9(x) the excitations. In addition, the abbreviations tp(x)-\-(p(— x) = 
2<f>(x), and 9(x) - 9(0) = Q(x) are used. 



2.4 The phase sensitivity 

We will use the phase sensitivity, i. e. the reaction of the system to a change in the 
boundary condition, to determine the localized-delocalized transition numerically for 
systems with finite size. We model the different boundary condition via a magnetic 
flux, penetrating the ring of the spinless fermions. The effect of the magnetic flux 
results in an additional phase of the hopping terms, tj — > \tj\e tei , with — tt < <dj < n. 
The energy levels depend only on the total flux, $ = arg (lljL^j). In particular, 
we determine below the energy difference between periodic (cjv = Co, $ = 0) and 
anti-periodic (cn = —cq, ^ — tt) boundary conditions, AE = (— ) Nf (E(0) — E(n)). 
Here the factor (— l) Nf cancels the odd-even effects resulting from the change in the 
ground-state for odd versus even particle numbers. We recall that, for a clean system, 
the ground-state energy and the finite-size corrections can be obtained from the Bcthc 
Ansatz Ipij . At half filling (and for odd particle number), the result, in the Luttinger 
regime, is given by |2j| 

E N ^)-Ne 00 = -^(^-3g^y (18) 

where is the energy density in the thermodynamic limit. Thus NAE = irvg/2, 
independent of N, for the metallic state. In an insulator, on the other hand, we find 
localized levels and the system cannot react to a twist in the boundary condition, i. e. 
NAE is expected to decrease with system size. Considering the above model, in which 
the contributions of the local potential and the weakening of the neighbouring bonds 
are taken into account separately, the localized levels do not split off the continuous 
spectrum p(|, and the phase sensitivity is an appropriate observable. 
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Fig. 1 Phase sensitivity as a function of interaction for a system with one site-impurity. 
The dashed-dotted lines correspond to Eq. (|20j), and the dashed lines to (|2l|). 



3 Comparison of site- and bond-impurities 

3. 1 Local potential 

A local magnetic field, which corresponds to a local potential, is the well-known and 
most studied example of a site-impurity. In this case, the free motion of the fermions 
inside the ring is restricted by the back-scattering at the impurity (±kp — * ^fkp). As 
discussed by Kane and Fisher the impurity becomes transparent for an attractive, 
and completely reflective for a repulsive interaction, according to the renormalization 
group equation 



de 
dl 



(19) 



In other words, the impurity strength scales either to zero or to infinity, for g > 1 and 
g < 1 , respectively. The fact that an impurity scales to zero for repulsive interaction is 
confirmed by the scaling of a weak link p2^ . For strong eo the effective hopping is given 
by t t = 4i 2 /e . The final result for the phase sensitivity is obtained by perturbation 



theory with respect to the defect strength, see 1 21 and [[27j , and is given by 



NAE = 



irvg 



N 
iVo 



1-9 



(20) 



in the weak scattering limit, and by 
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in the strong impurity limit, where Nq « 2 is a cut-off corresponding to a momentum 
cut-off of the order of the Fermi momentum. As a summary, we show numerical data 
in Fig. [l]. The most important result is that an intermediate defect strength scales 
to zero for attractive interaction, and to infinity for repulsive interaction, rather than 
scaling to an additional intermediate fixed point. 



3.2 One modified link 

We begin by checking whether the above expressions for the weak and strong potential, 



Eqs. ( [20D and (|2l|), can be applied to the case of a modified link, as should be expected 
from the treatment in [§2j or (2^]. We therefore calculate the phase sensitivity for a 
weakly modified link, St = 0.2, and a weak link, S t — 0.8. The results are shown in 
Fig. 0. However, the above expressions do not apply: Instead of Eq. ( |20| ) we find for 
the weakly modified link 

/ at \ (i— a)/2 

™9 or I N 



NAE=^-25 t {-) . (22) 

This is actually the scaling according to Eq. (|l4|). Therefore, Eq. ( p"2| ) does not 
appropriately describe this impurity type. Similarly, for the weak link, we do not 
recover Eq. (pl|), but instead 

*A2* = 4(1 -*)(—) . (23) 

In this case N — > N/2, as we connect two chains of length N/2. 

Next we verify that the generic link-impurity can be written as a sum of a hopping- 
and an interaction- impurity, as the bosonized Hamiltonians, Eqs. ([l4]) and (|l^), sug- 
gest. As shown on the 1. h. s. of Fig. ||, where the three cases: S t = S and Sy = 0, 
St = Sy, and 5y = § and St = 0, are compared, this conjecture is indeed confirmed. 
The qualitative behaviour is the same for all defect types. By comparing two different 
system sizes, we can fix the transition point to V — in all cases. 

If the impurity is so strong that the sign of the hopping or the interaction changes, 
S > 1, interaction dependent maxima occur, in addition, for the generic link- impurity, 
see the r. h. s. of Fig. |3[ 



3.3 Bon d-imp uri ty 

For the non-interacting system (V — 0) with one bond-impurity, it is easily seen that 
the back-scattering contributions, i. c. the Fourier components with q = n in the half- 
filled case, cancel each other as is the case for two site-impurities at an odd distance. 
For this reason, we presume that the phase sensitivity increases slightly with system 
size according to NAE oc (1 — 2bN /N) also in the interacting system. For repulsive 
interaction the results of Eggert & Affleck || show that this kind of perturbation is ir- 
relevant and the system remains metallic. For an attractive interaction, the behaviour 
is expected to be similar. Accordingly, the results for integrable, triangle like defects 
show that the phase sensitivity is slightly reduced for non-interacting fermions in 
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Fig. 2 Phase sensitivity as a function of interaction for a system with a hopping impurity. 
The straight line corresponds to the clean system, Eq. (|l8|). The dashed lines correspond to 
Eq. @, and the dashed-dotted lines to Eq. @. 
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Fig. 3 Left: Phase sensitivity as a function of interaction for a system with one site-impurity 
of type t, tV, and V. Right: Details of the phase sensitivity as a function of interaction for 
a system with one strong link-impurity. The straight line corresponds to the clean case, the 
other lines are connecting the data points. 
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Fig. 4 Left: Phase sensitivity as a function of interaction for a system with one bond 
impurity of varying strength, the system size is N = 40. The lines correspond to NAE(b = 
0)(1 — 2blo/N). Right: Phase sensitivity as a function of interaction for a system with one 
bond impurity of strength b = 0.5 (N = 40: +, N = 60: X, N = 80: A), and b = 0.9 
(N = 40: □, N = 60: □, N = 80: □). The lines correspond to fits according to Eq. @. 

the presence of such (transparent) impurities. In particular, a repulsive interaction 
leads to stronger suppression of the phase sensitivity than an attractive one. 

The numerical data presented in the left part of Fig. [I] confirm this assertion of 
an only slight reduction: clearly the system behaves metallic for small b. However, by 
increasing the impurity strength, strong deviations are visible, already for b = 0.3 at 
V = —1.9. The value of the interaction, where this strong decrease is seen, moves to a 
larger interaction when increasing the impurity strength, until the maximum reaches 
the repulsive regime for b > 0.7. Using the scaling applicable for the weak link case, 
Eq. (pl|), t t — > tf S = ttN 1 ' 1 / 9 , a scaling relation for the two-weak-link case, given by 
the Hamiltonian in Eq. (|i~7|), can be derived, 

(1 - b) -> (1 - 6) eff oc (1 - 6)iV 1 -9/ 2 - 1 /(4s)^ (24) 

According to this scaling relation, the weak links increase for V > —1.2. Thus, for 
an attractive interaction, there is an interaction value for which the phase sensitivity 
becomes independent of system size, while it scales to zero when further lowering 
the interaction down towards V — — 2 (where the first-order transition takes place 
in the clean system). Nevertheless, the system remains metallic in this region, but 
with a strongly reduced Drude weight. This behaviour, see the r. h. s. of Fig. [|, is 
different from the behaviour of an integrable defect (where a strong reduction arises 
for a repulsive interaction). 
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Fig. 5 Phase sensitivity versus interaction for a clean and a distorted dimerized system; the 
system size is N = 48. The lines correspond to Eq. (Eq). 



4 Combination of dimerization and impurity 

4.1 A barrier in a dimerized system 

We combine now two perturbations and study first a site-impurity in the dimerized 
system. For a first impression, we calculate the energy levels and the ground-state 
energy of of non-interacting particles with an alternating hopping, following the treat- 
ment in j28|, but adding one potential scatterer of strength e. The ground-state energy 
is given by 

k " I ^2 

E = -2t ^ Y cos2 ka + T^2 + " 2 sin2 ka K E ( K ) ( 25 ) 

with E(k) the second complete elliptic integral, where k is given by k 2 = (1 — u 2 )/(\ + 
e 2 (N)) ~ 1 — u 2 — e 2 (^V); e(N) — e/2N. A single impurity in the non-interacting 
alternating chain was also discussed in [p9[ . Thus the phase sensitivity is given by, 
compare | p0| : 

/ f N V~ 9 

NAE = NAE(0,0) - 2N ^u 2 + e 2 (N) - . (26) 

Figure shows numerical data in comparison to this formula, with No w 4 according 
to the cell doubling in the dimerized state. 

The length dependence presently is not as clear as for a single impurity or for 
the clean dimerized system: For g < 1 the system is completely localized, the phase 
sensitivity decreases with system size, while for g > 2, the system appears to be de- 
localized. In between, for 1 < g < 2, however, a characteristic system size N 2 = 
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Fig. 6 Phase sensitivity as a function of interaction for a system with u — 0.1 and e = 10. 
For V > —1 the phase sensitivity decreases with system size, and for V < —1.4 it increases. 
In between, it increases from TV = 24 to TV = 48, but decreases from TV = 48 to TV = 100. 



e 2 (g — l)/[4u 2 (2 — g)] can be defined. For TV smaller than N c , the phase sensitivity 
increases with TV similar to the behaviour in the delocalized phase. Increasing the sys- 
tem size above TV C , however, the phase sensitivity decreases, indicative of localization. 
This intermediate regime is shown in Fig. ^, where we plot the phase sensitivity for a 
set of parameters which clearly shows the described behaviour. 

An impurity with weak or intermediate strength enlarges the delocalized region 
of the dimerized model, as already shown in |3(i[| . For a strong barrier and a strong 
dimerization, a new feature is observed: In addition to the enlargement of the de- 
localized region a peak in the phase sensitivity is found in the localized regime, for 
e > 3 and u > 0.1. The characteristic behaviour is shown in Fig. ^, demonstrating 
that the increase for a value of Vpeak ~ — 1 is already by one order of magnitude for 
TV = 100. The assumption that the impurity becomes irrelevant at V pca k is consistent 



uy = 








uv = u t 
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Vpcak 


V e 
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Vpcak 


V c 


0.1 


10 


(-1) 


-1.3 


0.1 


10 




-0.9 


0.15 


10 


-l.i 


-1.4 


0.2 


10 




-1.3 


0.2 


10 


-1.2 


-1.6 


0.3 


10 


-1.0 


-1.45 


0.25 


10 


-1.2 


-1.7 


0.5 


10 


-1.05 


-1.65 


0.3 


10 


-1.2 


-1.8 


0.7 


10 


-1.1 


-1.75 



Table 1 Values of V^eak and V c for uv = 0, and uv = u t . The latter case is included for 
completeness, compare the discussion in subsection 2.2, even though not mentioned separately 
in the text. 
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Fig. 7 Typical plot of the phase sensitivity versus interaction for a system with strong 
dimerization and a strong impurity. 



with the initial increase of the phase sensitivity at this value. Since the effect of the 
dimerization is weakened by an impurity, the following (by further lowering the inter- 
action down to the value V c ) decrease of the phase sensitivity can be explained by the 
enhancement of the effective dimerization, caused by the diminished influence of the 
impurity. Lowering the interaction further, the dimerization becomes irrelevant, too, 
and the system delocalizes. Table |j] shows the characteristic interaction values for this 
model. 



4.2 Dimerization and bond-impurity 

While we analyzed above a single site-impurity in a dimerized ring, we now concentrate 



wo 



on a single bond-impurity plus dimerization. Based on the results of Sec. 3.3 
expect that for a weak distortion the system still shows the phase transition at V — 
—a/2, with the delocalized region extending to V = —2. Increasing the strength 
of the bond-impurity, we also expect the unusual behaviour - as described above - 
for strong attractive interaction to occur in the dimerized system as well, because 
in this interaction regime the dimerization is irrelevant. The first conjecture can be 
confirmed numerically for strong dimerization but small b. The numerical data for 
stronger distortion, i. e. for increasing b, show a more complex behaviour of the 
phase sensitivity, see Fig. |^. The first observation is, similar to a site-impurity, that 
the transition to the metallic state occurs at weaker interaction strength, e. g. for 
u = 0.03 already at V = —1 (compared to V ~ —1.5 in the undistorted case). Also, 
the second conjecture is supported by the numerical data: For an interaction larger 
than V w —1.5 the influence of the bond- impurity drives the systems again to a 
metallic state with a strongly reduced phase sensitivity. In between, —1.5 < V < — 1, 
the system flows to the free metallic fixed point when increasing the system size. 
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Fig. 8 Phase sensitivity as a function of interaction for a dimerized (u = 0.03) system with 
a strong bond-impurity, b = 0.7; the system sizes are N = 20 and N = 60, respectively. For 
comparison, the data of the clean dimerized system, again for N = 20 and 60, are included. 



5 Summary 

Our numerical studies of a one-dimensional spin chain, which is equivalent to a system 
of spinless interacting fermions, with various distortions show that both types of site- 
impurities are irrelevant for an attractive interaction (using the fermion picture), and 
relevant for a repulsive one. The description within a first-order perturbation treat- 
ment can be applied in the case of the site-impurity. The bond-impurity is irrelevant 
for repulsive interaction, and leads to a strong suppression of quantum coherence - but 
not to a localized ground-state for attractive interaction near the first-order transition. 
The scaling could be determined using the bosonization technique. In the dimerized 
system, the numerical results show that dimerization and any kind of impurity weaken 
each other. A single barrier in a dimerized chain leads to an enhancement of the dclo- 
calized region compared to the clean dimerized system. In addition, a sharp peak (as 
a function of the interaction parameter) , increasing the phase sensitivity by one to two 
orders of magnitude for small systems, N w 100, is found for strong perturbations. 

A bond-impurity in the dimerized system is irrelevant for a small impurity strength. 
A stronger impurity leads to a complicated phase diagram, especially for an attrac- 
tive interaction. The combination of dimerization and impurity therefore leads in all 
cases to a shift of the localization-delocalization transition of the dimerized system 
to a weaker interaction strength. By tuning the values of dimerization and impu- 
rity strength, different metallic windows in the insulating system can be opened. At 
an interaction of about V w —0.5, . . . , — 1 a bond-impurity turns the system into the 
metal. For a strong attractive interaction, V « —1.5, . . . , —2, where the bond-impurity 
reduces the Drude weight drastically a strong site-impurity leads to derealization, es- 



C. Schuster and U. Eckern, Quantum spin chains with defects 



15 



pecially in the strongly dimerized system. 
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